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Abstract. A recently proposed phase-estimation protocol that is based on measuring 
the parity of a two-mode squeezed vacuum state at the output of a Mach-Zehnder 
interferometer shows that sub-Hcisenberg sensitivity is obtained (PM Anisimov et 
al, PRL 104, 103602 2010). This protocol works in the case of infinite number of 
parity measurements. Here we consider the use of photon-number-resolving detectors 
to measure parity and we use Bayesian analysis to characterize the performance of 
the phase estimation in such a scheme but with finite number of measurements. We 
have found that phase estimation has non-zero bias at or 7r/2 phase. Yet there is an 
in-between region where the bias is negligible. Our phase estimation scheme beats the 
Heisenbcrg limit and saturates the Cramer-Rao bound. 
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1. Introduction 



Phase estimation is the primary tool for optical metrology. There are several 
experimental ways to estimate phase. Coherent light based interferometry is most 
commonly used but its sensitivity for phase estimation is limited by the shot-noise (SN) 
PQ. This is not a problem in the case of limitless resources or in the case of samples 
that can withstand large doses of radiation. However, this is a problem otherwise, and 
one has to resort to interferometry with quantum states of light, such as NOON states 
[2], and measuring parity [31 H] in order to achieve sub-shot-noise or even Heisenberg 
limited (HL) sensitivity of phase estimation. 

A recently proposed phase estimation scheme beats the HL in the case of infinite 
number of parity measurements [6]. This scheme measures the parity of the output 
signal in a Mach-Zehnder interferometer (MZI), Fig.HJ with two-mode squeezed vacuum 
(TMSV) input. It turned out that this particular scheme using TMSV input has sub- 
Heisenberg sensitivity even with linear phase evolution. This is due to the fact that the 
photon number uncertainty is greater than the mean photon number inside the MZI. 

Squeezed vacuum generation in an optical parametric amplifier (OPA) of up to 12 
dB of quadrature squeezing has been achieved experimentally [7]. This corresponds to 
a mean photon number in both modes of the TMSV of up to seven photons. Hence, 
the parity, even or odd photon number, of the state can be measured with existing 
photon-number-resolving detectors such as transition edge sensors [8j [9]. However, a 
photon-number-resolving detector does not provide a mean value of the parity signal 
after a single measurement, which means that a phase measuring experiment must be 
repeated multiple times. 

Our paper applies Bayesian analysis to this phase-estimation scheme when photon- 
number-resolving detectors are used to infer the parity signal. We use the parity signal 
at the output of the MZI to estimate the unknown phase 9. Our analysis shows that, 
although phase estimation is biased near the phase origin and at ir/2, there is an in- 
between interval where unbiased phase estimation is possible. In this interval, phase 
sensitivity saturates the Cramer- Rao bound (CRB), beats the HL, and remains below 
the SN limit. Unlike the recently proposed phase estimation scheme, we achieve these 
results with a finite number of parity measurements. 



2. Model 



Parity based phase estimation was originally introduced to quantum optics by Gerry in 
Ref. [ID] . The parity signal (II) for our phase estimation scheme, 



\J\ + n(n + 2) sin 2 6>' 

was obtained in Ref. [6] . It depends on the unknown phase 9 inside of the MZI and the 
mean photon n in the TMSV state used. Thus, knowing n and (II), the unknown phase 
9 can be determined. 
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Figure 1. A schematic view of a Mach-Zchnder interferometer (MZI). Two-mode 
squeezed vacuum states are generated at the input of the MZI by an optical parametric 
amplifier. We measure the parity signal at the output of the MZI with a photon number 
resolving detector such as a transition edge sensor. 
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Figure 2. Convergence of the parity signal (blue) to its expectation value of (II) = 0.93 
(red) in a finite number of M = 10 4 numerical runs. A different sample of this 
convergence will evolve differently due to the probabilistic nature of the scheme, but it 
will always converge to its expectation value. In this case the photon number is fi = 3 
and the unknown phase is 9 = 0.1. 



Each parity measurement returns either an even or odd outcome with probabilities 
P e and P . Thus, there is uncertainty in measuring the parity signal for our phase 
estimation scheme, so it requires averaging the outcomes of many parity measurements. 
The expectation value of a state's parity is (IT) = 1 ■ P e + (—1) ■ P Q . Since, P e + P a = 1, 
we can take advantage of the following expressions for the probabilities of an even and 
an odd outcome respectively: 

Pe = \ (l + (fi)) (2) 

Po = \ (l - (A)) (3) 

We use the probability of an even photon number, P e from Eq. (j2J), to numerically 
generate a value for the parity signal. Figure |2] shows how the parity signal (blue) 
converges to its expectation value of (ft) = 0.93 (red) in a set of M — 10 4 numerical 
runs for the case of an input photon number of n = 3 and an unknown phase 9 = 0.1. 
Repeating this procedure with the same parameters results in a slightly different parity 
evolution due to the probabilistic nature of the scheme. The uncertainty of the signal 
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Figure 3. The Bayesian probability density function (PDF) from Eq. (|5|) is evaluated 
with mean photon number n = 3 and numerical runs M = 200 for various number 
of outcomes to. We rescaled the heights for easier comparison of their widths and 
maxima. As the number of even outcomes decreases from to = M to m = 0.75M, the 
maximum of the PDF shifts towards <fi = ir/2 with a consequent broadening of the 
distribution. This broadening reduces the accuracy of estimating phase. 



leads to an uncertainty in the measurement of the unknown phase 6. We can improve 
the certainty of unknown phase estimation by repeating our measurement many times 
so that it converges to the expectation value. 

The Bayesian analysis of an even or odd parity detection at the output of the MZI 
provides a probability for the unknown phase 9 to be in the interval [</>, + dtp] given all 
prior observations. This analysis is based on the update rule for the probability density 
function (PDF) given by Bayes theorem: P(0|output) oc P (output |0)P pi . ior (0), so that 
each consecutive observation modifies the PDF and improves phase estimation. Hence, 
the update equation given by Bayes theorem for an even or odd outcome is: 



Particular outcomes in our numerical model are independent thus the Bayesian PDF 
after M runs with m even outcomes is: 



Figure [3] presents use of this update rule in the case of M = 200 runs with n = 3 input 
photons for different number of even outcomes: m = 150, 175, 200. The heights of each 
have been rescaled for the ease of comparison. 

One can see two effects that reduced number of even outcomes has on the PDF. 
The PDF becomes broader and shifts its maximum towards <fi = it/ 2. This results in 
higher phase uncertainty and thus reduction in the sensitivity of phase estimation near 
6 = 7r/2. This coincides with predictions of Ref. [6] that the best phase sensitivity is 
achieved in the vicinity of 9 = 0, where the parity of the output state is predominantly 
even. 

As M goes to infinity the PDF becomes 



P(0|{e,o})ocP({e,o}|0)P(0). 



(4) 




(5) 




(6) 
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Figure 4. Core phase dependence of Bayesian PDF for n — 3 in the limit of infinitely 
long averaging. A well defined maximum is present only near expected phase </> = 
and our unknown phase 9 — which explains good sensitivity at the vicinity of the 
phase origin. 
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Figure 5. Slices of the theoretical probability density obtained with Bayesian analysis 
over long averaging. We show various unknown phase's 9 with an input photon number 
of h — 3. Low values of 8 have well defined maximum, which provides a scheme for 
the fast and accurate estimation of the unknown phase 9. 

We use a phase that maximizes the PDF as our maximum likelihood estimator 
for an unknown phase 9. Success of this approach is governed by a core function 
P^ 9 ' n \<j),n)Po ^ e ' n \<f),n) that is presented in Fig. H] as a function of two phases, 9 and 
</>, for n — 3. Figure [5] shows slices of this core function. Only slices of the core function 
for low values of the unknown phase 9 demonstrate a well-defined maximum. Larger 
values of 9 yield ill defined maxima causing reduced accuracy and slow convergence 
of our scheme. A well defined maximum provides a scheme for the fast and accurate 
estimation of the unknown phase 9. 
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Figure 6. This distribution of phase estimations in 1000 trials shows uncertainty of 
phase estimation due to a finite number of numerical runs M for 9 = 0.08 and n = 3. 
Each phase estimation is obtained by inputting the number of even states m counted in 
M = 200 numerical runs into Eq. ([5]) and then finding its maximum. This distribution 
has the mean value of </> = 0.078 and a standard deviation of A0 = 0.022 which falls 
well within the unknown phase of 8 = 0.08. 



3. Numerical results 

Equipped with theoretical predictions, we now implement our phase estimation scheme 
numerically assuming that parity signal is measured with photon-number resolving 
detectors that detect even outcomes with the following probability: 



This comes directly from Eq. (j2J) which depends on an unknown phase 9 and input 
photon number n. In a single phase estimation, we repeat the process of outcome 
generation M times and determine a particular number of even outcomes. Using our 
Bayesian update rule from Eq. (J5J) with m even outcomes, we calculate the probability 
distribution for an unknown phase 9. Our maximum likelihood estimator concludes that 
the unknown phase is equal to the phase at the maximum of this probability distribution. 

Normally, phase sensitivity is estimated by measuring the width of a Bayesian 
PDF. In our case however, the Bayesian PDF is asymmetric and has several maxima in 
a 2-7T interval. Therefore, we analyse the statistics of observed maxima in the interval 
<ft G [0,7r/2] for a phase estimation instead. To determine uncertainty of our phase 
estimation, we repeat phase estimation 1000 times with M = 200 numerical runs in 
each. Figure [6] shows the distribution of different phase estimations in the case of 
unknown phase being 9 = 0.08 and n — 3. The distribution of phase estimations have 
a mean value of 4> = 0.078 with a standard deviation of A<p = 0.022 which falls well 
within the expected phase of 9 = 0.08. 

As we ran our phase estimator with input intensity n = 3 for several unknown 
phases in the interval 9 G [0,7r/2], we found the presence of bias, which is statistical 
favouritism that causes misleading results, defined as Bias = (ft — 9. 




(7) 
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Figure 7. Bias of phase estimation for n = 3 as a function of the unknown phase 9 
for increasing number of numerical runs M and 1000 trials in each phase estimation. 
Bias is negligible in the interval 9 6 [0.08, 0.6], and thus provides us with an accurate 
phase estimation scheme. 



Figure [7] demonstrates high bias in estimating the unknown phase near the 
boundaries of the interval. This is shown for a various number of numerical runs M 
which demonstrates how more measurements improve our scheme. For an unknown 
phase 9 G [0,0.08], our phase estimation exhibits a large bias toward zero. Likewise, 
when the unknown phase approaches n/2 our scheme again becomes biased. This leaves 
an optimal phase interval of 9 € [0.08, 0.6] that offers unbiased phase estimation. 

The interval for unbiased phase estimation depends on the input intensity n. For a 
larger n = 7, our phase estimator shows that the minimum value of the unbiased phase 
interval remains relatively unchanged; however, the length of the interval reduces due 
to signal localization near the phase origin. Namely, the right boundary of the interval 
reduces from 9 = 0.6 to 9 = 0.3. 

4. Phase sensitivity 

The uncertainty of phase estimation A0 quantifies the phase sensitivity of the scheme. 
In the case of conventional phase estimation with coherent laser light and intensity 
difference measurement, phase sensitivity is shot-noise limited, A<p > 1/vS. This 
sensitivity improves with increasing input intensity as well as the number of numerical 
runs. However, ultimate phase sensitivity is Heisenberg limited, > 1/ y/ Mn 2 , that 
has the same dependence on the number of numerical runs M but faster dependence on 
input intensity. 

Phase estimation with TMSV and parity detection is capable of beating the HL 
due to a greater photon number variance in the state. Phase sensitivity for this scheme 
is limited from below by the CRB [6]: 



A0 




(8) 



^Mn{2 + n) cos 9' 
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Figure 8. The standard deviation of phase estimation with an input photon number 
of n — 3 for increasing number of numerical runs M, We found that more phase 
estimates reduces the standard deviation. As 9 increases towards 7r/2 the standard 
deviation becomes large. This explains why bias is high when 9 > 0.6. 



Hence, an optimum phase sensitivity, A0 = 1/ ^ Mn{2 + n), obtained near 9 = is 
sub-Heisenberg. 

Our presented approach of measuring parity with photon-number-resolving 
detectors is however biased near 9 = 0. Thus, Bayesian phase estimation with photon- 
number resolving detectors can only be obtained in the low bias interval of 9 G [0.08, 0.6]. 
Figure M shows the standard deviation, A<p, of our phase estimation scheme. Increasing 
the number of runs in a phase estimate improves the phase sensitivity when 9 is near 
the origin. This is because phase uncertainty scales as A0 = c/yM, where the 
proportionality constant c depends on the unknown phase 9 and input photon number 
n. As we calculate the unknown phase 9 in various regions of bias with n = 3, we find 
that the standard deviation is low near the origin. 

As an example, consider 9 = 0.01 which is in a region of high bias. As we increase 
the number of numerical runs in a phase estimate, we find that the standard deviation 
remains constant and cannot be fitted to a c/yM function in order to determine the 
proportionality coefficient c (see Fig. [9]). The same behavior is observed in the other 
high-bias interval closer to 7r/2. 

In the following, we focus on unknown phases in the interval suitable for unbiased 
phase estimation with parity measurement using photon number resolving detectors. 
Our data does have the expected l/\/M dependence. Hence, we can compare the 
proportionality coefficient obtained from the fitting procedure with a value expected 
from the CRB as well as from the SN limit. 

Figure [10] presents the standard deviation of phase estimation as a function of 
the number of phase estimates for an unknown phase of 9 = 0.3 and n = 3. In this 
case, we find that the proportionality coefficient is ctmsv — 0.69. As we compare this 
proportionality coefficient with the SN limit, csn = 0.58, and the CRB, ccrb = 0.62, 
obtained from Eq. (jSJ), we see that the sensitivity is in close proximity to the CRB (red). 
With an unknown phase of 9 = 0.08 located closer to the phase origin in our unbiased 
interval we can improve this result. For a given unknown phase, the sensitivity of 
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Figure 9. For a single phase estimate we see how the standard deviation changes 
for increasing numberical runs M with n = 3 and unknown phase 9 = 0.01. We 
compare our results with that of the Cramer-Rao bound (CRB) and shot-noise (SN) 
limits. Since 9 = 0.01 is in a high biased region, the data does not have the l/VM 
dependence. We can obtain this dependence with the use of an unknown phase in the 
unbiased region as shown in Fig. 1101 
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Figure 10. The standard deviation versus the number of numerical runs M in the 
case of unknown phase 9 = 0.3 and photon number n = 3. It exhibits the expected 
1/y/M scaling with a proportionality coefficient, ctmsv = 0.69. It does not beat either 
ccrb = 0.62 or csn = = 0.58, but we can use an unknown phase of 9 = 0.08 

which is at the lower end of our unbiased interval to significantly improve this result. 



phase estimation could be improved by varying the mean photon number in the input 
TMSV state. This is analysed for 9 = 0.3 and the results are presented in Fig. [TT| 
where the phase sensitivity is characterized by ctmsv which removes the dependence on 
M and focuses on input intensity dependence. One can see that our previous choice 
of n = 3 is an optimum input intensity for this value of the unknown phase. For 
lower input intensities, phase sensitivity is in close proximity to the CRB. As the input 
intensity increases, the standard deviation of the current scheme becomes much worse 
than predicted by the CRB. Therefore, we investigate lower unknown phases in our 
unbiased interval. The best sensitivity of our phase estimation protocol is expected near 
the phase origin. Performing phase estimation in the unbiased interval offers significant 
improvement for phase sensitivity. With an unknown phase of 9 = 0.08 shown in 
Fig. H2J one can see that phase sensitivity of our scheme is nearly optimal and beats 
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Figure 11. The constant of proportionality, ctmsVj versus the number of photons in 
the input state for an unknown phase of 9 = 0.3. The case of n = 3 has the smallest 
standard deviation for phase estimation based on the scheme presented here (blue). 
The demonstrated sensitivity is in close proximity to the CRB in red. The SN limit in 
black and the Heisenberg limit (HL) in purple are present for comparison. 
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Figure 12. The constant of proportionality, ctmsVj versus the number of photons in 
the input state for an unknown phase of 9 = 0.08. Our numerical experiment (blue) 
does much better than SN (Black), almost as well as the CRB (red), and does as well 
as the HL (purple). 



the SN limit for an extended range of input intensities. Thus, parity measurement with 
photon-number-resolving detectors is close to the limiting performance for our scheme 
and could be used to demonstrate sub-Heisenberg limited phase estimation. 

5. Conclusion 

Phase estimation protocols that rely on parity measurement benefit from photon 
number-resolving-detectors that can be used to infer parity of the state. Here, we 
considered a particular phase estimation protocol that is based on TMSV input and 
parity detection at the output of the MZI. Use of photon-number-resolving detectors 
means that measurements must be repeated multiple times with Bayesian analysis 
applied to all outcomes. 

Our scheme shows that phase sensitivity saturates the Cramer-Rao bound and 
beats the Heisenberg limit with the use of a finite number of experimental runs. We 
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discovered that our maximum likelihood estimator is biased near the phase origin where 
the best sensitivity is expected. As a result, the standard deviation of the estimated 
values does not reduce with increasing number of phase estimations but stays constant. 
However, small values of the unknown phase in our unbiased interval 9 G [0.08, 0.6] allow 
for unbiased phase estimation. Phase sensitivity behaves as expected and remains in 
close proximity to the Cramer-Rao bound. Consequently, the phase sensitivity remains 
sub-shot-noise limited for a broad range of input intensities as long as phase estimation 
is performed near the origin. 

Acknowledgments 

We would like to acknowledge the Intelligence Advanced Research Projects Activity and 
the National Science Foundation for support. 

6. References 

[1] Jonathan P. Dowling. Quantum optical metrology the lowdown on high-nOOn states. Contemp. 

Phys., 49(2):125-143, 2008. 
[2] Hwang Lee, Pieter Kok, and Jonathan P. Dowling. A quantum rosetta stone for interferometry. 

J. Mod. Opt, 49(14-15):2325-2338, 2002. 
[3] Christopher C. Gerry and Jihane Mimih. The parity operator in quantum optical metrology. 

Contemporary Physics, 51(6):497-511, 2010. 
[4] R. A. Campos, Christopher C. Gerry, and A. Benmoussa. Optical interferometry at the heisenberg 

limit with twin fock states and parity measurements. Phys. Rev. A, 68:023810, Aug 2003. 
[5] C.C. Gerry and P.L. Knight. Introductory quantum optics. Cambridge University Press, 2005. 
[6] Petr M. Anisimov, Gretchen M. Raterman, Aravind Chiruvelli, William N. Plick, Sean D. Huver, 

Hwang Lee, and Jonathan P. Dowling. Quantum metrology with two-mode squeezed vacuum: 

Parity detection beats the heisenberg limit. Phys. Rev. Lett., 104(10):103602, Mar 2010. 
[7] Moritz Mehmet, Henning Vahlbruch, Nico Lastzka, Karsten Danzmann, and Roman Schnabcl. 

Observation of squeezed states with strong photon-number oscillations. Phys. Rev. A, 81:013814, 

Jan 2010. 

[8] B Cabrera, RM Clarke, P Colling, AJ Miller, S Nam, and RW Romani. Detection of single 

infrared, optical, and ultraviolet photons using superconducting transition edge sensors. App. 

Phys. Lett, 73:735-737, 1998. 
[9] K.D. Irwin and C.C. Hilton. Transition-edge sensors. In Christian Enss, editor, Cryogenic Particle 

Detection, volume 99 of Topics in Applied Physics, pages 81-97. Springer Berlin / Heidelberg, 

2005. 

[10] Christopher C. Gerry. Heisenberg-limit interferometry with four-wave mixers operating in a 
nonlinear regime. Phys. Rev. A, 61:043811, Mar 2000. 



